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We develop an analytical procedure to compute all relevant physical properties of scalar field 
oscillons in models with quartic polynomial potentials: energy, radius, frequency, core-amplitude, 
and lifetime. We compare our predictions to numerical simulations of models with symmetric and 
asymmetric double well potentials in three spatial dimensions, obtaining excellent agreement. We 
also explain why oscillons have not been seen to decay in two spatial dimensions. 



INTRODUCTION 
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It is well known that a real scalar field under the influ- 
ence of certain nonlinear potentials supports long-lived, 
spatially-localized, time-dependent configurations known 
as oscillons [T] |2]. Even though there are no topologi- 
cal or nontopological conserved charges associated with 
such objects, they exhibit remarkably long lifetimes, far 
exceeding naive estimates from linear approximations to 
the theory [5]. This property has prompted much in- 
terest in studying oscillons in more detail in two [H |S] 
and three [51 spatial dimensions. Structures with similar 
qualitative behavior have been observed in many differ- 
ent physical systems, from vibrating grains [^ to acoustic 
oscillations in the sun [H]. 

In the past two years, it has been shown that oscil- 
lons are not limited to real scalar field models; they have 
been found in 2d Abelian-Higgs models [9 , and in the 
Standard Model [TD]. They may also play a role in cos- 
mology [11], possibly having an impact on the dynamics 
of infiationary reheating, in symmetry-breaking transi- 
tions, and in false vacuum decay |12j . In Id de Sitter 
models, oscillons not only form from near-thermal initial 
conditions, but survive to retain approximately 50% of 
the total energy of the universe [13] . 

So far, most oscillon properties, such as their energy, 
radius, and lifetime, have been obtained numerically. It 
has been shown that they have a near-constant energy 
and a near-periodic motion with frequency just below 
the mass of vacuum excitations [51 [T3]. In spite of this 
progress, it is fair to say that a proper analytical under- 
standing of oscillons is still lacking. 

In the present work, we will derive all physically rele- 
vant properties of scalar oscillons in quartic polynomial 
potentials. Starting from a simple ansatzior their profile, 
we develop a nonperturbative analytical procedure to ob- 
tain their energy, lifetime, and minimum and maximum 
values of their frequency, radius, and core amplitude. We 
compare our analytical results to numerical simulations 
obtaining excellent agreement. 



We begin with the Lagrangian for a d-dimensional, 
spherically-symmetric real scalar field in fiat spacetime. 
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where a dot denotes a time derivative, and Cd — 
27r''/2/r((i/2). Inspired by the numerical solutions, which 
show that an oscillon is well approximated by a Gaus- 
sian undergoing near-periodic motion in its amplitude 
and with an effective radius oscillating with small am- 
plitude about a mean value [2j [3], we approximate the 
oscillon solution as 



cj){t, r) = A{t) exp 
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where A{f) = <j){r = 0, t)—(j>y is the displacement from the 
vacuum, (/),„. This approximation works particularly well 
in d = 2 and d — 3, although it clearly fails to reproduce 
the oscillon's large r behavior, where (j){r) ~ exp[—mr]. 

Since oscillons have been observed in a diversity of 
quartic polynomial potentials, we write V{4>) as [TS] 
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where the gj's are constants and V{(f>y) is the vacuum 
energy. Substituting eqs. |2] and [3] in eq. [T] 
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with 



n— 3 



d/2 



V{A) = Ua' + Y^(1) V"(0.M", (5) 
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where V"{<j>y) = d"V{cl)y)/df' and lj^ = V"{<j>^) + d/R'^ 
is the linear frequency. (Primes will also denote par- 
tial derivatives with respect to (/).) For potentials with 
V"{(f>v) > 0, V{A) is the potential energy of a nonlinear 
oscillator with at least one equilibrium point (at A = 0) 



for any value of R. For conservative motion we can write 
the energy as, 
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where ^max is the positive turning point (Amax > 0) of 
A{t). 

We also write the expression for the frequency lo as- 
sociated with the conservative motion between the two 
turning points, 
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where T is the period of oscillation, A ~ 
[2E/cR~2V{A)y^, CR = (^/2)'*/2i?'i, and A„,i„ is 

related to A^ax by ^(Amin) = ^(^max)- 

It was noted in [3J [13] that the location of A^ax rel- 
ative to the inflection points of V^(^) is an accurate in- 
dicator of oscillon existence. This can be understood 
by noting that since the inflection points separate re- 
gions of opposite curvature in V{A), the particular region 
probed by ^max indicates the degree to which the field 
configuration experiences the stabilizing nonlinearities of 
the potential V{(t>). Following [15], we therefore adopt 
"(^"(Aniax) as an indicator of oscillon stability: increas- 
ingly negative (positive) values lead to increased stabil- 
ity(instabihty). Defining /(Amax,i?) = V^"(^max), 
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In figure O] we show the (parabolic) level curve corre- 
sponding to /(^,nax,^) = 0. Since values inside (out- 
side) the parabola have / < (/ > 0), it is easy to show 
that, for a given R, there exists a stable region provided 
that R > Rmin where [TS] 
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As the figure indicates, when R = i?min, there is only 
one value of A^ax located in the stable region, given by 
Ao = -(4/3)^/2-^"7F^^. 



LIFE STORY OF AN OSCILLON 

In figure [T] the continuous and dashed lines show 
curves of constant energy (eq. pi as a function of R and 
^max for a symmetric double-well potential. 



Vi<t>) = 02 _ ^3 ^ ^4/4^ 
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that is, 51 = 0, 32 = 2, gs = -6, 54 = 6, (/)„ = 
in eq. ^ [Quantities were made dimensionless with 




Maximum Amplitude 

FIG. 1: Two curves of constant energy, Eosc ~ 41.3 (continu- 
ous line), and -Battract = 37.7 (dashed line), together with the 
parabolic level curve J(y4max,-R) = 0. Vertical dashed line 
locates the asymptote of / = 0. 



r'M 



{v9'^'Pv)~^ and the primes have been sup- 
pressed.] We also plot the parabola I{R, A^ax) = 
marking the boundary between the stable region (/ < 0, 
above the curve) and the unstable region (/ > 0, below 
the curve) . Note that the presence of an oscillon in the 
unstable region does not necessarily mean that it will im- 
mediately decay; the degree of stability is given roughly 
by the vertical distance of a point to the curve / = 0. 

Using the curve / = in fig. [T] we can describe the tra- 
jectory modeling the onset of an oscillon from an initial 
configuration starting at a somewhat arbitrary point A. 
On its way to an oscillon, the configuration will radiate 
excess energy until reaching point B. To determine the 
location of point B, note that the onset of the oscillon 
phase is marked by a minimum radius. (In what follows 
it will be clear that this must be true if radiation radi- 
cally slows down.) Taking the minimum radius given by 
eq. [9] and using that the oscillon stage should begin in 
the stable region, fixes the location of point B to A — Ao 
and R = Rmin, as shown in figure [l] As we shall see, this 
is confirmed by numerical results. For the potential of 
= 2.42 and Ab ^ Aq = (4/3)^/2 = 1.54, in 
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excellent agreement with the inset in figure [2] 

Since our model approximates an oscillon as having 
constant energy, we input these values into eq. [6] to 
obtain the oscillon energy, -Eosc^ 
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FIG. 2: The maximum amplitude and radius as a function of 
time for a configuration with Ro = 2.86 and Ao = 2 in the 
symmetric double well of eq. 10 (continuous lines), and for 



2.16.3^1^ 



with 



an asymmetric double well [V{(j)) — | 

Ao — 2 and Ro = 4 (dashed lines). The analytical predictions 
at the points {Ab, Rb) and {Ad, Rd) are indicated by arrows 
and dots. The insets show the minimum radius and frequency 
for the symmetric double well. 



where the numerical value given is for d — i and the 
potential of eq. 10 Compare with numerical results [2 
[3], where, roughly, 42 < E^sc < 45 during the oscillon's 
lifetime. 



a quartic potential. 
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where As was obtained above and the numerical result 
is for the potential of eq. 10 in d = 3. 

To obtain the final radius Rd , we evaluate eq. [TT] at 
oint D: Er,.r ^ {t: I2Y/'^R^V{Ad,Rd). Solving for 



point D: 
Rd, 
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where giAo) = V"Al/2 + J2l^^{2/n)'^/^V'''Al/n\. In 
d = 3, one gets a cubic equation for Ro- For the po- 
tential of eq. [To] with E^sc — 41.3 and A]j = 0.84, we 
obtain, Rjy — 3.43. In figure [2] we compare numeri- 
cal and analytical values for {Ab,Rb) (left arrows and 
dots) and {Ad,Rd) (right dots) for this potential (con- 
tinuous lines) and for an asymmetric potential (dashed 
lines). [The numerical simulation measures R by fitting 
a Gaussian to the oscillon and reading off the radius.] 

Since we know Amax(i) throughout the oscillon's life, 
we can use the potential V{A) to compute ^min and then 
use eq. [t] to obtain uj{t). For example, using ^max = 
A]j gives the decay (or critical) frequency, ujd = Wcrit = 
1.38 < Wniass = ^"(01.) = V2- In the inset of figure 
l2]we compare analytical (dot) and numerical values. It 
should be noted that this model correctly predicts that 
the oscillon frequency tends to a maximum (as observed 
in, for example, |14] ) as it approaches decay. 



In addition to the above estimate for ii^osc, we can de- 
rive an absolute lower bound for the oscillon energy. This 
is found by noting that the dashed line in fig. [T] shows 
a critical energy, denoted -^attract, below which the en- 
ergy curves do not probe the stable region. Using the 
/ = condition, one can eliminate R from the expression 
of the energy (eq. pi) , which becomes a function only of 
^max- Remarkably, this function has a minimum, the 
asymptote energy of figure [T] For the potential of eq. [TO] 

^attract = 37.7 {d = 3) and ^attract = 4.44 (d - 2). 

Having arrived at point B (after radiating excess initial 
energy) the oscillon proceeds leftward along its line of 
constant energy, ii^osc (Traveling rightward would cause R 
to tend to zero, which is unphysical). When the oscillon 
passes point C, its distance from / = begins to increase 
and its instability grows. To estimate the point of oscillon 
decay we note that, at point D, the vertical distance to 
the / = curve becomes infinite. Mathematically, point 
D is found by taking the limit i? — > oo in the expression 
7 = and solving for Au, the decay amplitude. For a 
cubic potential this gives. Ad = -{3/2)^^^{V" /V"). For 



OSCILLON LIFETIME 

Due to the remarkable near-periodicity of an oscillon, 
its spectral function in frequency space displays a series 
of narrow peaks, with the dominant one just below the 
mass frequency 5 . Only a small radiating tail penetrates 
the region with lo > Wmass- As in ref. [5 , we model 
the power spectrum of the oscillon core amplitude (p(0, t) 
around its dominant frequency peak with a Breit-Wigner 
profile. Since the power is proportional to the square of 
the amplitude of a given frequency component, we write 



a^{uj) = K 



2 



(14) 



where 0^(10) is the square of the amplitude of the fre- 
quency component of 0(0, t) with frequency lj, and K is 
a constant to be determined. Given that the amplitude of 
this primary frequency peak is approximately that of the 
oscillon itself (because the first peak is narrow and holds 
most of the power) a^((jJosc) = K f» ^max; the oscillon's 



core amplitude. Since the constant F setting the width of 
the resonance peak is related to the (inverse) timescale 
associated with the instantaneous radiation rate of the 
oscillon, it is natural to let F = ^, where E is the oscil- 
lon energy. 

Most of the (small-amplitude) radiation leaking from 
the oscillon is contained in a mode with wavelength com- 
mensurate with the size of the oscillon - that is, with 
wavelength '^ 4i?. The frequency of this wave is there- 
fore Wrad = ^ujI^, + {2tt/AR)\ 

Now consider an infinitesimally thin shell of radius R 
around the oscillon with volume CdR'^^^dr, which is filled 
with the outgoing radiation wave (w = Wrad) traveling 
with speed Wrad- Within this thin shell, the wave can be 
approximated as a Id plane wave of constant amplitude; 
hence its energy density is ^rad'^rad/^; where ^rad is the 
amplitude of the radiation wave in the shell. Since ^^ad 
is given by a^^oji-nd), the total radiation energy contained 
in the thin shell is a^(a;rad)^^^Crf-R''^^(ir. The radiation 
travels outward at a speed Urad = Wrad/fcrad, and dr = 
Vrnddt. The amount of energy lost by the oscillon per 
unit time, — _B, is therefore 



Finally, note that eq. 



E = A^-uj^^^CdR''-\,-.d 



1 dE 
\ 2 dE 
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If the oscillon is long lived, Wiad — ^osc S> dE/dE. In this 
case, and writing 4i = ^-s, we have 

' ^ dE dt E 



& + AaE'^ = 0, 



(16) 



where a = 2(1 - WoscMad)^[^^Cd-R'* ^Urad] ^■ 

To integrate eq. [16] we take a to be constant, a good 
approximation for an oscillon. Then, 
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where E^o is the asymptotic energy as i — > oo and Ei the 
initial energy. The oscillon decays at energy E^, > E^o 
in a time riifc given by 
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If Ei '^ E]j , the lifetime will be approximately zero, as ex- 
pected. When Ei ^ Ejj, the lifetime becomes indepen- 
dent of Ei and tends to a maximum, r,„axiifo — ^ ^ ^^ . 
To evaluate this expression, note that Ejy is simply the 
plateau energy obtained in eq. |11[ For Eqo, being the 
minimum energy the oscillon may possess, it is most nat- 
ural to choose Eaa — ^/attract- For the other parameters 
in a, we analytically calculate the average of each quan- 
tity over the stable phase of the oscillon's life {C -^ D 
in figure nj). For the potential V{(j}) = cjp — cj)^ + \(j)^, we 
obtain Tmaxiifc ~ 6300, which is quite accurate. 



'a{E — -Eattract) • 



16 admits the solution, E = 

Thus, HE ^ i^attract, S -> 0. In 

d = 2, one can show numerically that the lowest energy 
oscillon has E^sc — -^'attract = 4.44. We thus see why 
Tinaxiito — > cxD in d = 2: the oscillon reaches its attractor 
value and stops radiating. 

We have analytically computed all relevant physical 
properties of scalar field oscillons in quartic polynomial 
potentials, confirming that oscillons are indeed attractors 
in field configuration space and that their lifetimes can be 
estimated by studying their main radiating mode. Our 
results suggest why d — 2 oscillons seem to be stable: 
they reach their attractor state before decaying. We in- 
tend to extend our approach to oscillons in models with 
gauge fields and with more complicated potentials. This 
work was partially supported by a National Science Foun- 
dation grant PHY-0653341. 
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